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Abstract 

We propose a protocol for creating entanglement within a dissipative circuit QED network architecture that consists 
of two electromagnetic circuits (cavities) and two superconducting qubits. The system interacts with a quantum 
environment, giving rise to decoherence and dissipation. We discuss the preparation of two separate entangled cavity- 
qubit states via Landau-Zener sweeps, after which the cavities interact via a tunable "quantum switch" which is 
realized with an ancilla qubit. Moreover, we discuss the decay of the resulting entangled two-cavity state due to the 
influence of the environment, where we focus on the entanglement decay. 
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1. Introduction 

The coupling between an atom and the quantized 
electromagnetic field represents a paradigmatic setup in 
quantum optics that allows one to study the interaction 
of light and matter in a fully quantum way. In most 
related experiments, an atom is placed inside a cavity, 
from where stems the labelling "Cavity Quantum Elec- 
trodynamics" UJ. 

In the last years the implementation of such sys- 
tems in the solid state realm has been demonstrated 
13, S] . Depending on the experimental implementation, 
the atom is realized artificially with a superconducting 
circuit or a quantum dot, while a quantum LC-circuit as 
part of a transmission line or a photonic crystal, respec- 
tively, serves as cavity. This not only provides a widely 
enhanced tunability of parameters in an experiment, but 
also allows one to achieve the strong coupling limit be- 
tween the artificial atom and the cavity, where the inter- 
action exceeds the decay rates of the atom and the cav- 
ity by far. It has recently been shown that the coupling 
strength may even be enhanced by several orders of 
magnitude compared to conventional cavity QED exper- 
iments 101. Furthermore, circuit QED systems may be 
used for many experimental validations of quantum me- 
chanics, such as the quantum-non-demolition-like read- 
out of a qubit state |5], the generation of Fock states fl, 
the observation of Berry phases |7], or multiphoton res- 
onances 10]. Generation of entanglement between the 
qubit and the cavity via a Landau Zener sweep has been 
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theoretically proposed as well ITol] . From a funda- 
mental point of view these setups, realizing "quantum 
optics on a chip", give the opportunity to rediscover and 
improve the vast goals of quantum optics and to test ex- 
perimentally its concepts in the solid state. 

However, the mentioned experiments involve "only" 
a single qubit and one cavity. With a view to quantum 
computing, i.e. in order to realize a quantum register, 
it is highly desirable to couple several of these qubit- 
cavity systems to obtain a circuit QED network. This 
would not only allow for performing quantum opera- 
tions such as quantum gates, but also enable the cre- 
ation and distribution of entanglement throughout the 
network. Indeed, a setup to couple two cavities via an 



ancilla qubit has been recently proposed in 111 IL Il2ll . 
while in recent experiments jlSl Il4ll . two qubits have 
been entangled via one cavity. 



In this paper, referring to the setup of Ref. Ill 111 we 
suggest the transfer of entanglement between two qubit- 
cavity systems via an ancilla qubit, which requires en- 
tangling different cavities. We restrict ourselves to the 
case where the qubit-cavity entanglement is created via 
a Landau-Zener sweep; cf. Ref. llil lToll . 

Due to the interaction with a dissipative environment, 
entanglement within an open quantum system is typi- 
cally subject to decay via spontaneous emission. For 
the non-trivial class of states discussed here, this can 
even occur during finite time. This phenomena of "sud- 
den death of entanglement" has been investigated theo- 
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retically for a conventional cavity QED setup lIlSll and 
observed in an experiment as well lll6|] . 

The paper is organized as follows. In section |2] we 
shortly review the elementary description of a super- 
conducting qubit coupled to an electromagnetic quan- 
tum circuit. Further, we give a brief review on the cre- 
ation of entanglement between qubit and cavity by a 
Landau-Zener sweep. In section [3] we discuss the in- 
fluence of a dissipative environment on the final en- 
tangled qubit-cavity state. In section H] we consider a 
network architecture consisting of two qubit-cavity sys- 
tems dynamically coupled via an ancilla qubit, and dis- 
cuss how this dynamical interaction yields an entangled 
two-cavity state. Finally, in section |5] we investigate 
the finite-time disentanglement of that state under en- 
vironmental influence, depending on the details of state 
preparation. 

2. Qubit-cavity entanglement via LZ sweep 

Among many possible setups for circuit QED, one 
concrete realization is given by a Cooper pair box (CPB) 
that couples capacitively to a LC circuit, which acts as a 
harmonic oscillator |4]. The CPB is formed of a super- 
conducting island connected to a superconducting reser- 
voir by a dc SQUID. The effective Josephson energy 
Ej = cos(7r<l>/(I)o) can be tuned via an external flux 
O penetrating the SQUID, where Oq denotes the flux 
quantum. We assume that O is switchable within a suf- 
ficiently large time interval such that Ej-hvt, v >0. The 
capacitive energy of the CPB, Ef.\=jEc{N - Ng)^ is de- 
termined by the material-dependent charging energy E^, 
the number of Cooper pairs on the island, and the re- 
duced background charge Ng which is proportional to a 
voltage bias Vg. In the charging limit E^^ Ej and near 
the charge degeneracy point Ng=l/2, it is justified to de- 
scribe the CPB by its two lowest charge states |A^=0) and 
|A^=1). This allows the interpretation as a two-level sys- 
tem (TLS) or qubit for which the Hamiltonian pseudo- 
spin notation reads //qb= - j^hvtcr,. It possesses the 
eigenstates |T> = (|0> + 11»/ V2 and |i> = (|0> - 1 !»/ ^J2. 

Making use of the standard Rabi model the total 
qubit-oscillator Hamiltonian reads ||4J 

hvt 



-cr, + hgcTjtia' + a) + hQa'a 



(1) 



The second term refers to the coupling of the qubit to 
the fundamental mode of the transmission line, which is 
modeled as a harmonic oscillator with the usual bosonic 
creation and annihilation operators a' and a and en- 
ergy eigenstates |n), n=0, 1 , . . . , oo. Note that the usually 
performed Rotating Wave Approximation (RWA) in the 
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Figure 1: (Color online) (a) Adiabatic energy levels of the qubit- 
oscillator Hamiltonian (T) as a function of the Josephson energy which 
is swept at constant velocity such that £j = hvt. Here, g = 0.04CI 
and V = O.Oin^. (b) Dynamics for the probabilities \{l,\\i)'{t))\^ 
and Kt, 0|iA(0)P- The occupation of higher states is always less than 
0.01% (not shown). 



coupling Hamiltonian is not justified here since during 
almost the complete Landau-Zener sweep, the qubit and 
the harmonic oscillator are far detuned. 

If the effective Josephson energy is switched from a 
large negative to a large positive value, and assuming 
that both the TLS and the oscillator are initially in their 
ground state If, 0), the probability for the TLS to end up 
in the upper state can be computed exactly to read ifl^ 



(2) 



This clearly reminds one of the Landau-Zener formula 
for the bare two-level system 1 17 , 18 , 3, with the level 
splitting in the anticrossing at t-Q.lv (see Fig.[T]) now 
given by the coupling constant g. It should be noted 
that the probability Q does not provide any direct in- 
formation about the final oscillator state. However, due 
to the symmetry of the Hamiltonian ([T]i, every creation 
or annihilation of a photon is accompanied by a qubit 
flip, which restricts the resulting dynamics to the states 
It, 2«) and |J,, 2n+l). Furthermore, the "no-go-up" the- 
orem states that P-[,n^-\,m — for7«>« l(2(jl . This reduces 



the possible final qubit-oscillator states to the state 
|^(cx,))= VI - ^T-i IT, 0) + ^ c„U, 2« + 1) (3) 

with the normalization 2„ |c„|^ = 1. In the experimen- 
tally relevant limit of g <k Q, we numerically find that 
during the whole Landau-Zener transition, the state is 
well approximated by 



(4) 



with a(oo) -yjl - Pi^i , and /3{oo) - -y/^'t^i ^i, i.e. 
|ci I = 1 . We substantiate this by plotting a typical exam- 
ple for the characteristics of the LZ dynamics in Fig.[T] 
(b), where it becomes evident that finally only the states 
|t,0) and IX, 1) are significantly populated. This means 
a great simplification since it allows one to reduce the 
number of all relevant oscillator states to two. 

Then in particular, the qubit-oscillator state (|4]l is en- 
tangled and can therefore be mapped to the state of 
two entangled two-level systems. For future conve- 
nience it is useful to introduce the density matrix for the 
qubit-cavity system pq_c =|tA)(i/'l that reads, in the basis 

{|T,0>, U,0>, |T,l), 



Pq-c(t) 



|a(OP a*(t)/3(t) ^ 





aitWit) ^(f)|2 



(5) 



As a consequence, we can use the concurrence as a well- 
defined measure of entanglement for the two two-level 
systems. It is defined as C = max{/li - A2 - Aj, - A4, 0). 
The parameters A are the ordered eigenvalues of the ma- 
trix -y^(o-^, ® (ry)p(c7^y ® cr?) -y^ Here, the PauU 
matrices cTy and cr^, act on the qubit space and the re- 
duced oscillator space formed by the states |0) and |1), 
respectively. For the density matrix (|5]l, the concurrence 
is given by 



Cq_c = 2|a*(f)y6(f)l . 



(6) 



Thus, the amount of entanglement in the long-time limit 
is given by Cq-doo) - 2 ^J{1 - which just 

depends on the ratio Iv. As an example, we plot the 
typical dynamics for the concurrence in Fig. |2] from 
where it is evident that entanglement is created after 
passing the anticrossing. 

3. Influence of the dissipative environment 

In a realistic experimental scenario, the qubit- 
oscillator system has to be regarded as an open system. 



i.e. one that interacts with its environment, thus suffer- 
ing dissipation and decoherence. Dissipative effects in 
an electromagnetic circuit are characterized by a spec- 
tral density lio)) which, within a quantum mechanical 
description, can be modelled by coupling the circuit bi- 
linearly to the field modes of an electromagnetic envi- 
ronment [22 ]. For weak system-bath coupling, tracing 
out the bath degrees of freedom and following standard 
techniques, the bath can be eliminated within Bloch- 
Redfield theory 123^ ,24.1 yielding the quantum master 
equation 

h n n 



(7) 



with the anticommutator [A,B}+-AB + BA, the 
oscillator-bath coupling operator Q - a~+a and the 
operators Q-iYHs,, Q]/h and 



Q 



^ Jo Jo 



dojS(co)cos(a)T)Q(-T). (8) 



Here, S{a)) = I{aj) coth(ha)/2kBT) is the Fourier trans- 
form of the symmetrically-ordered equilibrium corre- 
lation function for a thermal environment at tempera- 
ture T. We assume the spectral density to be ohmic, 
i.e. I(a>) - cjy with an effective impedance y. The 
notation X(t) is a shorthand for the Heisenberg oper- 
ator u]^{t)XU{){t), where Uq is the system propagator. 
While the quantum master equation (jTjl is general, an 
explicit form with respect to the Hamiltonian ([T]i, to- 
gether with details about the numerical implementation 
can be found in Ref . Ii25i1 . 

3.1. Entanglement creation 

Due to the influence of the environment, in particular, 
the entanglement created between qubit and cavity is 
subject to decay during the Landau-Zener sweep. In the 
case of low temperatures, k^T- 0.01?iO, which is rele- 
vant for most experiments, our numerical investigations 
of the master equation ^ confirm the following sce- 
nario: Before reaching the first avoided crossing at time 
f < - Q/ V (see Fig. [TJ, the system remains in its ground 
state. Then at time t--Q.lv, it evolves into the super- 
position (IDi. Simultaneously, the bath causes coherence 
decay and, for low temperature, a relaxation from the 
state It, 1) to It, 0) occurs. We emphasize here that there 
is no transition |t, 0) |i, 0) since the qubit experiences 
an effective heat bath with a spectral de nsity sharply 
peaked at the oscillator frequency Q 
Thus, for large times t » Q/v, the spectral density at 
the qubit splitting hvt vanishes and, consequently, the 
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Figure 2: (Color online) Qubit-oscillator entanglement in terms of the 
concuiTence Cq-c for a Landau-Zener transition. Here, k^T=O.Olh£l, 
while g and v are as in Fig. \i\ For finite coupling strength y > (red 
line), the decay of the concurrence is approximated by expression 111! 
(dashed black line). For comparison, the concurrence is also given for 
the case of vanishing interaction with the environment (blue line). 



qubit is effectively decoupled from the bath. As a con- 
sequence, we numerically find that after the anticross- 
ing the state |i, 0) becomes populated in addition to the 
states |T,0) and \i, 1). Quantitatively, the corresponding 
density matrix takes the form [cf. Q] 



Pq-c(t) = 



a c \ 

z 



c* li ; 



(9) 



After the Landau-Zener transition, its entries can be 
well approximated by 



(10) 



where tlz = ^/v and Pi^i has been defined in Eq. (O. 
For the state (|9]l, the concurrence is now given by 

Cq_c = 2\c\ - 2 VPT-i(l-^T-i)e"''™- d D 

Hence the loss of entanglement is directly related to the 
loss of coherence between the qubit and the cavity and 
exhibits a simple exponential decay, see Fig. |2] 



4. Qubit-cavity networks 

In this section we discuss how entanglement may be 
created within a quantum network architecture consist- 
ing of several pairs of cavities and qubits. We shall 
make use of the setup proposed in Ref. 1.12.] together 



Figure 3: (Color online) Sketch of the quantum network under consid- 
eration. Qubits are represented by big dots and cavities by black lines. 
Qubit A (B) is only coupled to cavity A (B), whereas both cavities may 
effectively interact with each other by the ancilla qubit. 



with the state preparation via a Landau-Zener sweep 
discussed above, in order to entangle two cavities. We 
study the resulting dynamics in presence of dissipation 
and decoherence. The setup we have in mind is sketched 
in Fig. |3] and consists of two qubit-cavity systems. The 
cross geometry indicates that qubit A (B) is only cou- 
pled to cavity A {B) whereas the central qubit (ancilla) 
is coupled to both cavities and may be used as a switch, 
as we will explain below. The Hamiltonian of such a 
system can be approximated as 



H 



{ if) 



= 1,2 



-I- G(a'^ + a^)(a'g + a^) . 



(12) 



The first line of Eq. (fT2l i denotes the usual Rabi Hamil- 
tonians for the respectively coupled qubits and cavities. 
For simplicity, both cavities are supposed to have the 
same frequency Q. The second line describes the an- 
cilla qubit and its coupling to the cavities, whereas the 
last line accounts for the direct geometrical interaction 
between the cavities which will be referred to as "geo- 
metrical coupling" [11]. 

4.1. The entanglement protocol 

In order to entangle the cavities A and B, we pursue 
a "quantum switching protocol" in analogy to the pro- 
ceeding described in Ref. Ol ill : In order to reach the dis- 
persive coupling regime g.i/A <K 1 between ancilla and 
cavities with the detuning A - Wanc - the ancilla qubit 
is widely detuned from the cavity frequency 14J. On the 
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Figure 4: (Color online) Sketch of the quantum switch protocol, dis- 
playing the time-dependent qubit and cavity energies and the ancilla 
interaction: The qubits A and B undergo Landau-Zener sweeps with 
their respective cavities with velocities I'^i and vg. This leads to the en- 
tangled qubit-cavity state 0. Meanwhile, the cavities do not interact, 
i. e. gsw = (dashed line). From t = Tan to t = T^g, gsw > (solid 
line), leading to the entangled state )15t which is subject to decay for 
t > Toff. For further details, see text. 



Other hand, A<s;f2, which allows a rotating-wave ap- 
proximation with respect to the ancilla-oscillator cou- 
pling. Under this condition the ancilla gives rise to 
an effective "dynamical coupling" between the cavities 
which may be expressed as glcr,(a'^ag + a^a^)/A. The 
total interaction between both cavities including the ad- 
ditional geometrical coupling [last line of ( fT2b 1 is then 
given by 



^^int - gsw (a^ag + a^Gg) , gsw 



—cr, + G. 

A " 



(13) 



It is possible to "switch off" this interaction by setting 
gsw to zero. This requirement can be fulfilled varying 
A or g;t, or manipulating the state of the ancilla. The 
protocol is now performed as follows: At first stage, 
both cavities and qubits become entangled via Landau- 
Zener sweeps at different sweep velocities va and vb, 
while the ancilla is in the "off" state, gsw-0. After the 
Landau-Zener sweeps, the state of the total system is the 
product state 



Pcav,i(0 



(14) 



■ PqA-cA WPqB-cB, 

where Pq,A/B-cq,A/B are of the form At this point, 
there is no entanglement between any object of A and 
B. Once the state (fT4] i is achieved, the interaction is 
"turned on" towards a finite interaction strength gsw > 
at time t - Ton- This may be accomplished e.g. by chang- 
ing adiabatically the operational point of the ancilla, 
which modifies A |11]. In the following, correlations 
between the cavities A and B emerge by virtue of the in- 
teraction (fT3T l. At time t = Toff, the interaction is turned 
off again. As far as we are only interested in the en- 
tanglement dynamics of the cavities, we may trace out 



the qubit and ancilla degrees of freedom in (fT4] i. For 
Toff - Ton = '2-n/gsw the resulting density matrix pcav.f 
takes the form 



Pcav,f 



a 








5(Zl +Z2) 
j(Zl -Z2) 







•-(zi-zi) 



jiZl +Z2) 





^ 





is 



(15) 



in the basis {\QaOb}, |0a1b), IUOb), IUIb))- For simplic- 
ity, we have introduced the notations 

a - qaOb + oaZb + aaZA + ZaZb 
j6 = dAdB 

Z\ = {ZA+aA)dB 

zi - (zB + aB)dA- (16) 

The coefficients aA/B,dA/B and za/b are defined in 
Eq. ([Tol l, where the indices A and B refer to the re- 
spective cavities. This reveals the dependence of the 
final state pcav.f on the specific parameters of our partic- 
ular entanglement protocol. We point out that the most 
relevant parameters for the preparation of ( fTsT i are the 
Landau-Zener sweep velocities va and vb- The corre- 
sponding concurrence is given by 



Ccav.f = 2 max \ 0, ^ |zi 



■Z2\ 



(17) 



As a result, using the state preparation via Landau- 
Zener sweeps together with the quantum switching pro- 
tocol described above, we end up with both cavities be- 
ing in an entangled state that depends on the particu- 
lar parameters of the protocol. The individual entan- 
glement between the individual qubit-oscillator systems 
has thus been transferred via cavity-cavity entangle- 
ment. 

5. Finite-time disentanglement 

In the presence of a dissipative environment of each 
cavity, the entanglement between the two cavities is nat- 
urally subject to decay. Now the possibility of finite- 
time disentanglement in a quantum optical setup due to 
spontaneous emission has been reported in Refs. 1115 , 
bon for the special class of entangled two-cavity states 
such as pcav,f ; see Eq. ( fTsl l. In this case the concurrence 
C(f) decays within a finite time, depending on the ini- 
tial state and thus on the particular state preparation. 
This signifies a "sudden death of (non-local) entangle- 
ment" whereas the decay characteristics of local deco- 
herence is exponential. Finite-time disentanglement has 
also been observed in a recent experiment 1.16,1 . 




Figure 5: (Color online) Time evolution of the concun'ence Ccav(0 
starting with state Pcav.v Here, va = 0.0050.^, gA = gs = 0.04n 
and T = O.OOlfi. White color con'esponds to separable states. Here, 
I'"" = 0.0624n-. 



In the following we study the disentanglement of the 
state pcav,f ■ Again, as an important point, our proposed 
circuit QED setup and protocol would allow for this to 
be analyzed in the solid state realm, while Refs. llisllsoll 
refer to an optical cavity QED setup. 

We start with the effective Hamiltonian for the uncou- 
pled cavities 



^etf 
-"cav.f 



(18) 



Following section |2] we may restrict the description of 
the cavities to the two lowest Fock states. In order to in- 
vestigate the time evolution of the entangled two-cavity 
state pcav,f after the ancilla interaction has been set to 
zero at f = Toff, we solve the quantum master equation 



eff 



f,p]+{.lA{aA+ai)+-CB(aB+al))p{l9) 



where the Liouvillian X.a/b is defined as in Eq. (|7]l- 
Hence we assume both entangled cavities to interact in- 
dividually with vacuum noise, i.e. to interact with in- 
dividual environments. A scenario with two two-level 
systems with zero distance, i.e. coupled to the same en- 
vironment, was discussed in Ref. jsiil . 

In Fig. |5]we plot the concurrence Ccav(0 for the de- 
caying state pcav.f as a function of time and the sweep 
velocity v^, while keeping va constant. We find that the 
concurrence drops to zero at a finite time f > fd for a 
broad range of different values of vg. This indicates that 
the corresponding initial states become disentangled af- 
ter a finite time. However, this does not always hold: 
For vb > V™' we rather find exponential concurrence 
decay, i.e. = oo. Note that, in contrast to these results, 
the qubit-cavity entanglement always decays exponen- 
tially as discussed in Section[3] 



Our results comply perfectly with the disentangling 
behavior discussed by Yu and Eberly fisll . which un- 
derlines the generic nature of that scenario. 

6. Conclusions 

We have proposed a protocol for creating an entan- 
gled two-cavity state within a circuit QED network ar- 
chitecture. Thereby, we have considered a network of 
two superconducting qubits coupled each to an electro- 
magnetic quantum circuit. While there is no direct in- 
teraction between the qubits, the cavities are connected 
through a third "ancilla" qubit. Depending on its state 
and its detuning from the cavities, this coupling may 
be switched on and off. As first protocol step, one has 
to prepare two entangled qubit-cavity states, which we 
have suggested to be done via Landau-Zener sweeps. 
Note that there exist other possibilities for this purpose, 
for example by Rabi oscillations. Second, the interac- 
tion is "switched on" in order to entangle both cavities. 
At the same time we have considered the influence of a 
separate dissipative envkonment for each cavity. 

A major goal here is the preparation of an entan- 
gled two-cavity state such as (flST l. In the quantum op- 
tics community this class of states has been proposed 
to investigate finite-time disentanglement due to spon- 
taneous emission. Here we have shown that this can 
be done as well in the solid state realm, depending on 
the particulars of state preparation. The entries of the 
density matrix (flST l can be obtained by measuring cross- 
correlations between the cavities. A correspondin g cir - 
cuit QED experiment has been proposed recently ll32ll . 
Thus, we belief that our protocol will enable experimen- 
talists to study non-trivial dynamics of entanglement in 
quantum circuits. 
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